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Abstract 

Charmed (and bottom) hypernuclei are studied in the quark-meson coupling (QMC) 
model. This completes systematic studies of charmed (A^",E C ,S C ), and hypernuclei in 
the QMC model. Effects of the Pauli blocking due to the underlying quark structure of 
baryons, and the T, C N — A C N channel coupling are phenomenologically taken into account 
at the hadronic level in the same way as those included for strange hypernuclei. Our 
results suggest that the X+ + and hypernuclei are very unlikely to be formed, while 
the A+, H[l and hypernuclei are quite likely to be formed. For the hypernuclei, the 
formation probability is non-zero though small. A detailed analysis is also made about 
the phenomenologically introduced Pauli blocking and channel coupling effects for the S[? 
hypernuclei. 
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1 Introduction 



Based on the quark-meson coupling (QMC) model [1], we have recently initiated a systematic 
study of charmed and bottom hadrons in nuclear matter [2] . The results imply that the forma- 
tions of i?~-meson nuclear (atomic) bound states, and A+ and A& hypernuclei are quite likely. 
The formations of A+ and/or A& hypernuclei were first predicted in mid-1970s by Tyapkin [3], 
and Dover and Kahana [4]. Later on theoretical studies [5, 6] as well as possibility of experi- 
mental observation [7] were made. Further, we have studied [8] the A+ and A b hypernuclei by 
solving equations of motion for a system of hypernuclei, embedding a A+ or a A& into the closed 
shell nuclear core, within the Hartree, mean-field, approximation using the QMC model. The 
results again support the possible formation of A+ and A& hypernuclei. 

In this article we present the results for the S c and S c hypernuclei, as well as some of 
the A+ and A b hypernuclei studied recently [8]. This is an extension of the studies made for 
the strange hypernuclei [9], and a completion of our systematic study of charmed and bottom 
hypernuclei [2, 8] in the QMC model. The extension is straightforward, since the heavy baryons 
(as well as nucleon) are classified by the naive constituent quark model, and QMC can treat 
them on the same quark-based footing. 

One of the other important issues in this study is a role of the u and d light quarks 
in nuclear medium, in the charmed or bottom baryons containing u and/or d quarks. This 
concerns the partial restoration of chiral symmetry in nuclear medium for heavy baryons. 

Concerning the experimental possibilities, the approved construction of the Japan Hadron 
Facility (JHF) with a beam energy of 50 GeV, will produce charmed hadrons profusely and 
bottom hadrons in lesser numbers, but still with an intensity which is comparable to the present 
hyperon production rates. Although the study of possible existence of charmed and bottom 
hypernuclei may be in its infancy at present, it is clear that conditions for the experiments to 
search for such heavy baryon hypernuclei are now becoming realistic, and would be realized at 
JHF [10]. Thus, it is still meaningful to consider such heavy baryon hypernuclei at this stage 
for the future experimental studies. 

The QMC model (QMC-I), which is used in this study, has been successfully applied to 
many problems of nuclear physics, such as finite nuclei [11, 12], strange hypernuclei [9], nuclear 
matter [13], and other problems including hadronic properties in nuclear medium [14, 15, 16, 
17, 18, 19, 20]. (See also Refs. [21, 22] for earlier, different versions of QMC, and Ref. [23] 
for related work.) For example, the model was applied to the study of meson nuclear bound 
states [15, 16], J/ty dissociation in nuclear matter [18], D and D production in antiproton- 
nucleus collisions [19], and kaon properties in nuclear matter and kaon production in heavy 
ion collisions [20]. Furthermore, although only some limited studies for heavy mesons (not for 
heavy baryons) with charm in nuclear matter were made by the QCD sum rules (for J/\I/ [24, 25] 
and D(D) [26]), a study for heavy baryons containing charm or bottom quarks was made only 
recently [2] using the QMC model. Furthermore, recent measurements of polarization transfer 
in the 4 He(e, e'p) 3 H and 16 0(e, e'p) 15 N reactions performed at MAMI and Jefferson Lab [27], 
support the medium modification of the proton electromagnetic form factors predicted by the 
QMC model [17]. 

There is also another version of the QMC model (QMC-II), where masses of the meson 
fields are also subject to the medium modification in a self-consistent manner [28]. However, 
for a proper parameter set (set B) the typical results obtained in QMC-II are very similar to 
those of QMC-I. The difference is ~ 16 % for the largest case, but typically ~ 10 % or less. (For 
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the effective masses of the hyperons, the differences between the two versions of QMC turn out 
to be less than ~ 8 %. See also Fig. 2 in Ref. [15] for the differences in the scalar potentials for 
the 1] and uj mesons in 208 Pb nucleus.) Since we use QMC-I, the differences in these numbers 
may be regarded as uncertainties in the model. 

Certainly, the model has shortcomings that have to be improved eventually. Difficulties 
to handle the model will increase rapidly in handling the Hartree-Fock approximation even 
for nuclear matter [29], and further to include the Pauli blocking effect at the quark level. In 
addition, the T,N — AN channel coupling effect has not been implemented yet in a consistent 
manner with the underlying quark degrees of freedom [9] . It should be mentioned that in the 
case of larger mass number E hypernuclei, no narrow states have been observed experimentally, 
although |He was confirmed [30]. It is unlikely to find such (narrow) bound states in the present 
situation, where the experimental analysis [30] supports the suggestion made by Harada [31], 
that a large isospin dependent E-nucleus potential term exists and the 1/A (A:baryon number) 
dependence of that term reduces the likelihood of observing bound states for A > 5. This un- 
likeliness may also arise from the width of £ state due to strong X — A conversion. Furthermore, 
an application to double hypernuclei has not been attempted, although recently the existence 
was confirmed for small baryon number nuclei [32]. As an attempt to improve the model, the 
effect of short-range quark-quark correlations associated with the nucleon overlap was studied 
in Ref. [33]. With the addition of the correlations, the saturation curve for symmetric nuclear 
matter was greatly improved at high density. Despite the shortcomings of the model described 
above, there are some positive aspects of this model, its simplicity and successful applicability, 
that we feel confident that the QMC model will provide us with a valuable glimpse into the 
properties of charmed and bottom hypernuclei. 

In the QMC model, the interactions between nucleons are mediated by the exchange of 
scalar (<r) and vector {uj and p) fields self-consistently coupled directly to the quarks within 
those nucleons, but not to the nucleons as in Quantum Hadrodynamics (QHD) [34]. Thus, it 
can be systematically applied to study the properties of any hadrons in nuclear medium [2, 14] 
(if they contain u and/or d light quarks). We make use of this advantage of the QMC model, 
and study the charmed and bottom hypernuclei. We hope that the present study inspires some 
interest to measure the properties of heavy baryons in nuclear medium, at facilities like JHF, 
particularly at RHIC (high energy relativistic heavy ion facilities), and in very high energy 
experiments at CERN and Fermilab. 

The organization of this article is as follows. In Section 2, the relativistic formulation 
of the charmed and bottom hypernuclear system in the QMC model will be explained briefly. 
A mean-field Lagrangian density and equations of motion will be given in Section 2.1, while 
charmed and bottom hadrons in nuclear matter will be discussed in Section 2.2. In Section 3, 
the effects necessary for a realistic calculation for the charmed (bottom) hypernuclei will be 
described. Spin-orbit potential, and the effect of the Pauli blocking at the quark level will be 
explained in Sections 3.1 and 3.2, respectively. In addition, Y, C ^N — A Ct t,N channel coupling 
effect will be explained in Section 3.2. The results for the charmed and bottom hypernuclei will 
be presented in Section 4, and finally Section 5 will be devoted to summary and discussion. 
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2 Charmed and bottom hypernuclei in the QMC model 



In this Section, we briefly explain the mean-field equations of motion for a charmed or bottom 
hypernuclear system, and review some properties of the charmed and bottom baryons in nuclear 
matter. 



2.1 Mean- field equations of motion 



We consider static, approximately spherically symmetric charmed and bottom hypernuclei, i.e., 
closed shell nuclear core plus one charmed or one bottom baryon (heavy baryon) configuration, 
ignoring non-spherical effects due to the embedded heavy baryon. The existence of the heavy 
baryon inside or outside of the nuclear core (in particular, A& which interacts relatively strong 
with core nucleons) breaks spherical symmetry, and one should include this effect in a truly 
rigorous treatment. We have neglected this effect, since it is not expected to be so important 
for spectroscopic calculations [35, 36] for the middle and/or large baryon number charmed 
hypernuclei (but may not be true for A& hypernuclei), and beyond the scope of present study. 
However, we include the response of the nuclear core arising from the self-consistent calculation, 
which is a purely relativistic effect [36, 37]. Thus, we always specify the state of the heavy 
baryon in which the calculation is performed, because the core nucleus response is different due 
to the self-consistent calculation. (This self-consistent procedure will be applied also when the 
effective Pauli blocking and/or £ Cj b — A Ct b are introduced later.) 

We adopt the Hartree, mean-field, approximation. In this approximation the pNN tensor 
coupling contributes to a spin-orbit force for a nucleon bound in a static spherical nucleus, 
although in Hartree- Fock it also contributes a central force [11, 12]. (Furthermore, it gives no 
contribution for nuclear matter since the meson fields are independent of position and time.) 
Thus, we ignore the pNN tensor coupling in this study as in Ref. [9]. 

A relativistic Lagrangian density for a charmed or a bottom hypernucleus in the QMC 
model may be given by [9]: 
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c c 

'-QMC 



2 lv v 2 

J2 ^ C (f) [<7 • d - M*(a) - ( + gflfm + eQ c A(f) ) 7o ] ^ c (f), (1) 



where ipnif) ijPcif)) an d b(r) are respectively the nucleon (charmed or bottom baryon) and the 
p meson (the time component in the third direction of isospin) fields, while m a , m w and m p are 
the masses of the a, u and p meson fields. g u and g p are the u-N and p-N coupling constants 
which are respectively related to (u, (i)-quark-u;, g%, and (u, d)-quark-p, g q , coupling constants 
as g^ = 3g q and g p = g q [11, 12]. Note that these meson fields, a, uj and p, represent the 
quantum numbers and Lorentz structure which mediate the interactions among the nucleons 
in nuclear medium as in original QHD [34], corresponding, a <-> <f>o, uj <-> Vq and b <-> bo, and 
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they may not be directly connected with the physical particles, nor quark model states. Their 
masses in nuclear medium do not vary in the present treatment of QMC-I. Hereafter we will 
use notations for the quark flavors, q = u,d and Q = s, c, b. 

In an approximation where the a, u and p fields couple only to the u and d quarks, the 
coupling constants for the heavy baryon are obtained as g% = (n q /3)g U! , and = g p = g q p , with 
n q being the total number of valence u and d light quarks in the heavy baryon C . 1% and Qc 
are the third component of the heavy baryon isospin operator and its electric charge in units 
of the proton charge, e, respectively. The field dependent a-N and a-C coupling strengths 
predicted by the QMC model, g a (a) = g^ (a) and g%(a), related to the Lagrangian density of 
Eq. (1) are defined by 

M*(a) = M N -g a (a)a(f), (2) 
M*( a ) = Mc -gC(a)a(r), (3) 

where {Mc) is the free nucleon (heavy baryon) mass. Note that the dependence of these 
coupling strengths on the applied scalar field must be calculated self-consistently within the 
quark model [9, 11, 12]. Hence, unlike QHD [34], even though g^(a)/ g a (a) may be 2/3 or 1/3 
depending on the number of light quarks in the heavy baryon in free space, a = 0, (but only 
when their bag radii in free space are exactly the same), this will not necessarily be the case in 
nuclear matter. Explicit expressions for g^(a) and <7 CT (cr) will be given later in Eq. (14). From 
the Lagrangian density Eq. (1), a set of equations of motion for the heavy baryon hypernuclear 
system is obtained: 

[< 7 • d - M* N (a) - ( guW (i>) + g^btf + |(1 + r 3 N )A(f) frollMO = 0, (4) 

[i 7 • d - M*(a) - ( gZu(r) + g p I°b(r) + eQ c A(r) ) 7 o]^c(r) = 0, (5) 

/ _9 9x rdMtria), . . r dM%(a), r ,. 

(-V r 2 + ml)a(r) = -[ ^ > ]p s {f) - [ £ > \p c s {r), 

= g«C N (a)p s (r) + g c a C c {<j)p C s {?), (6) 
(-V r 2 + ml)uj{r) = gu) p B (r) + g%p%(r), (7) 

(-V r 2 + ml)b(f) = | P3 (r) + g c p if pg(fO, (8) 

(-V 2 r )A(r) = ep p (r) + eQ c p c B (r), (9) 

where, p s (f) (pf (f)), Pb{t) {Pb{^))i Pzif) an d Ppif) are the scalar, baryon, third component of 
isovector, and proton densities at the position r in the heavy baryon hypernucleus [9, 11, 12]. 
On the right hand side of Eq. (6), -[dM%{a) / da] = g a C N (a) and -[dM£(a)/da] = g°C c (a), 
where g a = g a (a = 0) and g% = g%(a = 0), are a new, and characteristic feature of QMC 
beyond QHD [34, 41, 42, 43]. The effective mass for the heavy baryon C is defined by 



n q g q a J ba d 3 x ij q (x)ij q (x) = -n q g q a S c (a) = [g° {?)<*] , (10) 



da a Jbag 

with the MIT bag model quantities, and the mass stability condition to be satisfied in obtaining 
the in-medium bag radius [1, 9, 11, 12]: 

M* c {*)= E ^^ + U R cfB, (11) 
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n;/2 + m;i%(n;-i) n2 , 

5c( } " n;(fi;-i) + m^/2' (12) 

^ = ^ + {R* c m$\ Q* Q = Jx% + (R* c m Q )*, m* q = m q - <#r(f), (13) 

C7 c (<t) = 5c(a)/5 c (0), = n,<&S c (0) = ^g„r c/N = ^g o S c (0)/S N (0), (14) 



dM* c 



dR c 



= 0. (15) 

Rc=R* c 



Quantities for the nucleoli are similarly obtained by replacing the indices, C — > N, in Eqs. (10) - 
(15). Approximating the constant, mean meson fields within the bag neglecting the Coulomb 
force, the wave functions, ipf(x)(f = u,d,s,c,b), satisfy the Dirac equations for the quarks in 
the baryon bag centered at a position r of the nucleus, (\x — f| < bag radius [15, 16, 20]): 

iT-^-K-^-^^Oi^W^^) = 0, (16) 

[«7 • d x -m Q }i) Q {x) = 0, (17) 

where, q = u or d, and Q = s, c or b. The constant, mean-field potentials within the bag 
are defined by V«{r) = g«a(r), V«{r) = &w{f) and V«{r) = ^fe(f), with g«, g« and ^ the 
corresponding quark-meson coupling constants. The eigenenergies, c u ^,q in units of 1/R* C are 
given by 

(^)=n; + ^(^(r)±^(r)), ^ = ^ (is) 

where, f2* and flq are given in Eq. (13). 

In the expressions of the MIT bag quantities Eqs. (11) and (13), zc, B, x q; Q, and m qj Q 
are the parameters for the sum of the cm. and gluon fluctuation effects, bag pressure, lowest 
eigenvalues for the quarks q or Q, respectively, and the corresponding current quark masses, zn 
and B (zc) are fixed by fitting the nucleon (heavy baryon) mass in free space. For the current 
quark masses we use (m Ujd , m s , m c , m b ) = (5, 250, 1300, 4200) MeV, where the values for m c and 
m b are the averaged values from Refs. [38] and [39], respectively, and these values were used 
in Refs. [2, 8, 16, 18, 19]. Then, we obtain the bag pressure B = (170 MeV) 4 , by choosing the 
bag radius for the nucleon in free space Rn = 0.8 fm. The quark-meson coupling constants, 
which are determined so as to reproduce the saturation properties of symmetric nuclear matter 
are, (g*,gl,g«) = (5.69,2.72,9.33), where g a = g% = 3<^(0) = 3 x 5.69 x 0.483 = 8.23 [12]. 
(See Eq. (14) with the replacement, C N.) The parameters z,-, and the bag radii Rj in 
free space (j = N, A, S, S, A+, S c , S c , A 6 ) and also some quantities calculated at normal nuclear 
mater density ps = 0.15/m~ 3 , are listed in Table 1, together with the free space masses [40]. 
(Nuclear mater limit will be discussed in Section 2.2.) 

At the hadronic level, the entire information on the quark dynamics is condensed into the 
effective couplings C N: c(cr) of Eq. (6). Furthermore, when Cjv,c(c) = 1; which corresponds to 
a structureless nucleon or a heavy baryon, the equations of motion given by Eqs. (4) - (9) can 
be identified with those derived from QHD [41, 42, 43], except for the terms arising from the 
tensor coupling and the non-linear scalar and/or vector field interactions introduced beyond 
the naive QHD. 

The explicit expressions for coupled differential equations to obtain various fields and 
baryon wave functions to describe a heavy baryon hypernuclear system, for the spherically 
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Table 1: The bag parameters, the baryon masses and the bag radii in free space [at normal 
nuclear matter density, p B = 0.15/m~ 3 ] Zj,Rj and Mj [R* and MJ] (j = N, A, A+, E c , S c , A 6 ), 
respectively. 



j 




Mj (MeV) 


Rj M 


M* (MeV) 


R* (fm) 


N 


3.295 


939.0 


0.800 


754.5 


0.786 


A 


3.131 


1115.7 


0.806 


992.7 


0.803 


£ 


2.810 


1193.1 


0.827 


1070.4 


0.824 




2.860 


1318.1 


0.820 


1256.7 


0.818 




1.766 


2284.9 


0.846 


2162.5 


0.843 




1.033 


2452.0 


0.885 


2330.2 


0.882 




1.564 


2469.1 


0.853 


2408.0 


0.851 


A, 


-0.643 


5624.0 


0.930 


5502.9 


0.928 



symmetric nucleus plus one heavy baryon configuration can be obtained by replacing the cor- 
responding indices from hyperon to heavy baryon (Y — > C) in Ref. [9]. 

2.2 Nuclear matter limit 

Here, we summarize the results and discussions given in Refs. [2, 9] for a limit of nuclear matter. 
New things in this subsection are some results given in Tables 1 and 2, and to inform the similar 
treatments (parameterizations) also work for charmed and bottom baryons as will be shown in 
Eqs. (21) and (22), which may not be so trivial. In the nuclear matter limit all meson fields 
become constant, and we denote the mean-values of the u and a fields by u; and a, respectively. 
Then, equations for the ZJ and self-consistency condition for the a are given by [1, 11, 12, 13] 

- = wrl i3k9(kF - k) ^3^f B ' (19) 

- - (20) 

where g„ = 3(^5^(0) (see Eq. (14) by the replacement, C — > N), Icf is the Fermi momentum, 
p B and p s are the baryon and scalar densities, respectively. Note that M^(a) in Eq. (20), 
must be calculated self-consistently by the MIT bag model, through Eqs. (10) - (16). This 
self-consistency equation for the a is the same as that in QHD, except that in the latter one 
has Cn(7j) = 1 [34]. Using the mean field value a, the corresponding quantity for the heavy 
baryon C, Cc(cx), can be also calculated using Eqs. (10) - (17) neglecting the effect of a single 
heavy baryon on the mean field value a, in infinite nuclear matter. 

It has been found that the function Cj(a)(J = N, A, S, S, A+, S c , S c , A 6 ) can be parametrized 
as a linear form in the a field, g a a, for practical calculations [9, 11, 12]: 

Cj(a) = 1 - a, x (g a a), (j = N, A, E, S, A+ S c , ~ c , A 6 ). (21) 

The values obtained for a,j are listed in Table 2. This parameterization works very well up to 
about three times of normal nuclear matter density ps — 3p (po — 0.15 fm -3 ). Then, the 
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effective masses for the baryon, j, in nuclear matter are well approximated by [9, f f , 12] 



n n 



1 " y (9*v)\ o, U = N, A, E, S, A+ E c , S c , A b ), 



(22) 



with n q being the number of light quarks in the baryon j. For the field strength g a a versus 
baryon density, one can find in Ref. [11]. 



Table 2: The slope parameters, aj (j = N, A, E, S, A+, E c , S c , A b ). 





xlO" 4 MeV" 1 


a 3 


xlO~ 4 MeV" 1 


UN 


8.8 




10.9 


a A 


9.3 


°A+ 


9.8 


a s 


9.5 




10.3 


a= 


9.4 




9.9 



Next, we discuss briefly the scalar and vector potentials of charmed and bottom baryons 
in nuclear matter [2]. The scalar (Vj) and vector (Vj) potentials for the baryon j, in nuclear 
matter are given by 

V> = m*-m v (23) 
Vi = n q V« + liV«, (24) 

where 1% is the third component of isospin projection of the baryon j. Thus, the vector potential 
for a heavy baryon with charm or bottom quark(s), is equal to that of the hyperon with the 
same light quark number in the QMC model. Calculated results for the scalar potentials for 
the baryons in symmetric nuclear matter are shown in Fig. 1, which is taken from Ref. [2] by 
omitting those for mesons. From the results it is confirmed that the scalar potential for the 
baryon j, Vj , follows a simple light quark number counting rule: 

V> c |vf, (25) 

where n q is the number of light quarks in the baryon j, and V S N is the scalar potential for 
the nucleon. (See Eq.(23).) It is interesting to note that, the baryons with charm and bot- 
tom quarks (S c is a quark flavor configuration, qsc), show features very similar to those of 
the corresponding strange hyperons with the same light quark numbers. Then, ignoring the 
Coulomb force although it is important in a realistic hypernucleus, we may expect that these 
heavy baryons with charm or bottom quarks, will also form some heavy baryon hypernuclei 
at this stage, as the strange hyperons do. (Recall that the repulsive, vector potentials are the 
same for the corresponding hyperons with the same light quark numbers.) 



3 Calculations 

In this section, we discuss parameters used in the calculations: spin-orbit force in the QMC 
model, the Pauli blocking effect at the quark level, and the effect of channel coupling, E Cj6 A^ — 
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A Cj j,iV. The parameters are the same as those used in the study of strange hypernuclei [9], except 
for those related to charmed and bottom hypernuclei given in Section 2, but the treatment is 
the same as that of Ref. [9]. More detailed discussions may be found in Ref. [9]. 

3.1 Spin-orbit force in the QMC model 

The origin of the spin orbit force for a composite nucleon moving through scalar and vector 
fields which vary with position is explained in detail in Ref. [11]. For the A, £ and S cases are 
explained in Ref. [9] . Although the spin-orbit splittings for the nucleon calculated in QMC are 
already somewhat smaller [11, 12], it was demonstrated that much smaller spin-orbit splittings 
were obtained for the A in QMC [9]. These will turn out to be even smaller for the heavy 
baryon hypernuclei. 

In order to include the spin-orbit potential approximately, e.g., for the A+, we add per- 
turbatively the correction due to the vector potential, — 2M J ^ r {^g^ u{rj^ I ■ s, to the single- 

particle energies obtained with the Dirac equation, in the same way as that added for the 
strange hypernuclei [9]. This may correspond to a correct spin-orbit force which is calculated 
by the underlying quark model [9, 11]: 



since the Dirac equation at the hadronic level solved in usual QHD-type models leads to: 



which has the opposite sign for the vector potential, g^ c oj{f). The correction to the spin-orbit 
force in the QMC model, appears due to the structure of the baryon (finite size). This may 
also be modeled at the hadronic level of the Dirac equation by adding a tensor interaction, 
motivated by the quark model [44, 45] . In addition, one boson exchange model with underlying 
approximate SU(3) symmetry in strong interactions, also leads to a weaker spin-orbit forces for 
the (strange) hyperon- nucleon (YN) than that for the nucleon-nucleon (NN) [46]. 

In practice, because of its heavy mass (M*+), contribution to the single-particle energies 
from the spin-orbit potential both with or without the inclusion of the correction term, turns out 
to be even smaller than that for the A hypernuclei, and further smaller for the A& hypernuclei [8] . 
Contribution from the spin-orbit potential with the correction term is typically of order 0.01 
MeV, and even the largest case is ~ 0.5 MeV, among all the heavy baryon hypernuclei calculated 
in this study. This can be understood when one considers the limit, M* + — > oo in Eq. (26), 
where the quantity inside the square brackets varies smoothly from an order of hundred MeV 
to zero near the surface of the hypernucleus, and the derivative with respect to r is finite. The 
results for the other heavy baryon hypernuclei are quite similar, and spin-orbit forces give only 
tiny contributions for the single-particle energies. Thus, the spin-orbit splittings are expected 
to be small. 

3.2 Effects of the Pauli blocking and channel coupling 

Below, we briefly explain the effects of the Pauli blocking due to the quark structure, and the 
channel coupling [9]. 




(26) 




(27) 
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Although the present treatment is based on the underlying quark structure of the nucleon 
and heavy baryons, it is missing an explicit inclusion of the Pauli blocking effect at the quark 
level among the u and d quarks in the core nucleons and the heavy baryons. For the strange 
hypernuclei, the effect of the Pauli blocking was included in a specific way at the hadronic 
level. We follow the treatment of Ref. [9], and include effective Pauli blocking and also the 
^c,bN ~ A-c,bN channel coupling effects, although the channel coupling effect is expected to be 
smaller than that for the EiV — AA, since the mass difference for the former is larger than the 
latter. Thus, the channel coupling effect that is included for the S c hypernuclei in the present 
study should be regarded as the limiting case. 

Following Ref. [9] , we assume that the Pauli blocking effect is simply proportional to the 
nucleon baryonic density (or u and d total light quark number density), although one could 
consider a more complicated density dependence. Then, the Dirac equation for the heavy 
baryon C, Eq. (5), may be modified by 

[i 7 • d - M* c {a) - ( \cPb{t) + 9^uj(r) + g p I°b{r) + eQ c A(r) ) 7 c#c(r) = 0, (28) 

where, Pb{t) is the baryonic density at the position f in the heavy baryon hypernucleus due 
to the core nucleons, and Ac is a constant, taken the same value to that determined for the 
strange hyperon. We note that the Pauli blocking effect associated with the light quarks 
in the heavy baryon should also lead to some repulsion for the nucleons. Contrary to the 
comment made in Ref. [9], the actual calculations in this study (and in Ref. [9]) include the 
effect of the modification made in Eq. (28) in a self-consistent manner. In practice, the above 
modified Eq. (28) has been inserted again into a system of coupled differential equations for the 
spherical static (heavy baryon) hypernuclear system, from Eq. (4) to Eq. (18), and calculations 
are performed again. Thus, the response of the core nucleons due to the modification in Eq. (28) 
for the heavy baryon, is included. 

For the A hypernuclei [9] the parameter corresponding to Ac in Eq. (28) was chosen to 
reproduce the empirical single particle energy for the lsi/2 in ^ 09 Pb, -27.0 MeV [48]. The fitted 
value for the constant was, Aa = 60.25 MeV (fm) 3 . Thus, we use the same value also for the 
A C) 6, i.e., Aa c6 = 60.25 MeV (fm) 3 . For the S c and H c hypernuclei, we can deduce the constants 
As c ,h c , corresponding to the effective Pauli blocking effect as, As c = Aa c , and A Sc = §Aa c , by 
counting the total number of u and d quarks in the heavy baryons. 

One may regard that this fitted value includes also the attractive A^N — > E^A channel 
coupling effect for the A Cj b single-particle energies, since the value fitted for the A hypernuclei 
includes the experimentally observed effect of the AiV — * £A channel coupling. As for the 
E C A — > A C N channel coupling effect, it must be included in addition to reproduce the relative 
repulsive energy shift in the single-particle energies for the S c hypernuclei, as was done for the 
£ hypernuclei. 

The £A — AA channel coupling was estimated using the Nijmegen potential in Ref. [9]. 
The effect for the X was included by assuming the same form as that for the effective Pauli 
blocking via AsPb(?")j and adjusted the parameter, As = Aa — > As ^ As, to reproduce this 
difference in the single-particle energy for the lsi/2 in g; 9 Pb, namely, -19.6 = -26.9(~ -27.0) + 
7.3 MeV. The value obtained for As in this way was, As = 110.6 MeV (fm) 3 . Thus, we use the 
same value, As c = 110.6 MeV (fm) 3 , although the effect is expected to be smaller considering 
the mass differences, ms — m^ ~ 77 MeV and ms c — m A + ~ 166 MeV [40]. As for the EN — AA 
channel coupling the effect was estimated to be negligible in Ref. [9], and we also neglect the 
effect of channel coupling on E c in this study. 
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4 Results for heavy baryon hypernuclei 



Before presenting the results, we summarize in Table 3 the parameters at the hadronic level 
fixed by the saturation properties of symmetric nuclear matter. Concerning the parameters for 
the a field, we note that the properties of nuclear matter only fix the ratio, (g a /m a ), with a 
chosen value, m a = 550 MeV. Keeping this ratio to be a constant, the value m a = 418 MeV 
for finite nuclei (hypernuclei) is obtained by fitting the r.m.s. charge radius of 40 Ca to the 
experimental value, r ch ( 40 Ca) = 3.48 fm [12]. (See also Section 2.1 for the value of g a in nuclear 
matter.) 



Table 3: Parameters at the hadronic level for finite nuclei [12]. 



field 


mass (MeV) 


3 2 /4tt (e 2 /47r) 


a 


418 


3.12 


00 


783 


5.31 


P 


770 


6.93 


A 





1/137.036 



First, in Tables 4 and 5, we summarize the calculated single-particle energies for 1 O, 41 Ca, 
49 Ca, 91 Zr and 209 Pb (j = A, A+, A 6 ) hypernuclei, together with the experimental 

data [47, 48] for the A hypernuclei. The results for the A+ and A b , and A hypernuclei are 
taken from Refs. [8], and [9], respectively for comparison. In the calculation, we have searched 
for the single-particle states up to the same highest state as that of the core neutrons in each 
hypernucleus, since the deeper levels are usually easier to observe in experiment. 

After a first look at the results shown in Tables 4 and 5 we notice: 

1. X+ + and 2+ hypernuclei are very unlikely to be formed in usual circumstances. Solid 
conclusions may not be drawn about !?++„+ Ca. Although results imply the formations 
of these hypernuclei, those numbers may be regarded within the uncertainties of the 
model and the approximations made in the calculation. Only a major alteration of the 
parameters may yield such bound hypernuclei. 

2. S° and hypernuclei may have some possibilities to be formed. However, due to a 
peculiar feature that the correct l«i/2 state is not found in ^>0, | 9 Ca and |oZr, one needs 
some caution to take this statement naively. (Discussions on this matter will be made 
later in detail.) A calculational modification could increase confidence in predicting such 
hypernuclei. 

3. A+,H° and A& hypernuclei are expected to be formed quite likely in a realistic situation. 
However, for the A& hypernuclei, it will be very difficult to achieve such a resolution to 
distinguish the states experimentally. The question of experimental observation is not 
tackled realistically. 

4. The Coulomb force plays a crucial role in forming (unforming) hypernuclei, by com- 
paring the single-particle energies among/between the members within each multiplet, 
(A,A+,A b ), (£°, £+,£++), and (~°, ~+) hypernuclei. 
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Table 4: Single-particle energies (in MeV) for ) 7 0, fCa and fCa (j = A, A+, S c , S c , A 6 ). 
Results for A+ and A b , and A hypernuclei are taken from Refs. [8], and [9], respectively. Ex- 
perimental data are taken from Ref. [47]. Spin-orbit splittings for the A hypernuclei are not 
well determined by the experiments. 



To 



~ Tr o 



17 O 

E++ U 



IP 



O (Expt.) 



17 D 



1P3/2 
lPl/2 



-12.5 
-2.5 (lp) 
-2.5 (lp) 



-14.1 

-5.1 

-5.0 



-12.8 

-7.3 

-7.3 



-10.7 
-10.2 



-7.5 
-4.0 
-3.6 



-7.9 
-3.5 
-3.5 



-2.1 



-19.6 
-16.5 
-16.5 





i°Ca (Expt.) 


^Ca 


A+ Ca 


jiCa 


£+Ca ^++Ca 


in C a 

'— 'c 




lSl/2 


-20.0 


-19.5 


-12.8 


-16.3 


-6.3 


-9.9 


-1.0 -23.0 


IPS/2 


-12.0 (lp) 


-12.3 


-9.2 


-13.2 


-4.1 


-6.7 


-20.9 


lPl/2 


-12.0 (lp) 


-12.3 


-9.1 


-12.9 


-3.8 


-6.7 


-20.9 


14/2 




-4.7 


-4.8 


-9.9 




-3.3 


-18.4 


2si/ 2 




-3.5 


-3.4 


-9.3 




-2.8 


-17.4 


14/2 




-4.6 


-4.8 


-9.4 




-3.3 


-18.4 











^t)Ca 




^ 9 ++ Ca 


-Jo Ccl 

— 'c 


1 9 +Ca 


A 9 (,Ca 


lSl/2 




-21.0 


-14.3 




-7.4 


-5.0 


-7.8 


-5.1 


-24.4 


1P3/2 




-13.9 


-10.6 


-7.2 


-5.1 


-3.4 


-4.1 


-2.8 


-22.2 


lPl/2 




-13.8 


-10.6 


-7.0 


-4.9 


-3.2 


-4.1 


-2.8 


-22.2 


14/2 




-6.5 


-6.5 


-4.0 


-2.3 




-0.9 




-19.5 


2si/ 2 




-5.4 


-5.3 


-4.9 






-1.5 




-18.8 


14/2 




-6.4 


-6.4 


-3.6 


-1.9 




-1.0 




-19.5 


1/7/2 






-2.0 












-16.8 
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Table 5: Single-particle energies (in MeV) for fZr and f & ¥h (j = A, A+, E c , S c , A b ). Results for 
A+ and A&, and A hypernuclei are taken from Refs. [8], and [9], respectively. Experimental data 
are taken from Ref . [48] . Spin-orbit splittings for the A hypernuclei are not well determined by 
the experiments. The values given inside the brackets for |oZr are those obtained by switching 

off both the Pauli blocking and the channel coupling effects simultaneously by setting As c = 0. 





i 9 Yb (Expt.) 


fZr 




91 7 r 

c 




91 7 r 


^oZr 

' 'c 


91 +Zr 


91 7 r 




-22.5 


-23.9 


-10.8 


"(-) 


-3.7 




-9.3 




-25.7 


1P3/2 


-16.0 (lp) 


-18.4 


-8.7 


-10.2 (-26.5) 


-2.3 





-6.6 


— 


-24.2 


lPl/2 


-16.0 (lp) 


-18.4 


-8.7 


-10.1 (-26.4) 


-2.1 





-6.7 


— 


-24.2 


14/2 


-9.0 (Id) 


-12.3 


-5.8 


-7.6 (-21.8) 


— 





-4.0 


— 


-22.4 


2si/ 2 




-10.8 


-3.9 


-8.1 (-23.0) 


— 





-3.9 


— 


-21.6 


14/2 


-9.0 (Id) 


-12.3 


-5.8 


-7.3 (-21.6) 


— 





-4.0 


— 


-22.4 


1/7/2 


-2.0 (1/) 


-5.9 


-2.4 


-5.1 (-17.1) 


— 





-1.3 


— 


-20.4 


2^3/2 




-4.2 


— 


-5.0 (-16.5) 


— 





-1.3 


— 


-19.5 


1/5/2 


-2.0 (1/) 


-5.8 


-2.4 


-4.7 (-16.8) 


— 





-1.4 


— 


-20.4 


2pi/ 2 




-4.1 


— 


-4.9 (-16.3) 


— 





-1.3 


— 


-19.5 


109/2 




— 


— 


-2.4 (-12.4) 


— 





— 


— 


-18.1 




l 08 Pb (Expt.) 


f>Pb 


f + 9 Pb 


|°o 9 Pb 


299 Pb 


^Pb 


|?, 9 Pb 

'c 


l° + 9 pb 




1 Si /O 

J -°l/2 


-27.0 


-27.0 


-5.2 


-7.5 






-6.7 




-27.4 


J-F3/2 




9^ 4 


4 1 


u.u 










9(S (S 


L Pl/2 


99 D fin) 


9^ A 


4 n 


yj.o 










9fi fi 


lu 5/2 


-17 (1 rfl 


-1 Q 1 


-9 4 


-5 3 






(J. u 




-95 4 


2 Si /o 




-17.6 










-3.3 




-24.7 


1 rf/o /o 

±u, 3/2 


-17 (Id) 


-19.1 


-2.4 


-5.1 






-4.0 




-25.4 


J 1/2 


-12 (1 f) 


-14.4 




-3.8 






-2.2 




-24.1 






-12.4 










-1.4 




-23.2 


1 A/9 
X J 5/2 


-12 (1 f) 


-14.3 




-3.5 






-2.3 




-24.1 


2j0i /? 




-12.4 










-1.5 




-23.2 


109/2 


-7.0 (1<?) 


-9.3 




-2.1 










-22.6 


107/2 


-7.0 (lg) 


-9.2 




-1.8 










-22.6 


1^11/2 




-3.9 














-21.0 


24/2 




-7.0 














-21.7 


24/2 




-7.0 














-21.7 


1^9/2 




-3.8 














-21.0 


3S1/2 




-6.1 














-21.3 


2/7/2 




-1.7 














-20.1 


3^3/2 




-1.0 














-19.6 


2/5/2 




-1.7 














-20.1 


3pi/ 2 




-1.0 














-19.6 


1*13/2 


















-19.3 
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We discuss below in detail the TP C hypernuclei single-particle energies, which have a pe- 
culiar feature. In some of the hypernuclei calculated, the lsi/2 state could not be found in 
the following sense. We have searched for the lsi/2 state by varying matching point to get the 
eigenenergy and the corresponding eigen wave function for the E° in steps 0.1 fm, from 0.2 
fm to 5.0 fm from the center of the hypernucleus. However, always a proper level with lower 
energy than those of the lpi/2, IP3/2 and 2s!/ 2 states could not be found, although a smaller 
matching point than that for the 2s\/2 state and the Dirac's k — —1 necessary to obtain the 
correct lsi/2 wave function were input. In such a case, the single-particle energy level became 
equal or higher than that of the 2si/ 2 state, which is also above the levels of the lpi/2 and IP3/2 
states. In the search procedure, we always checked that the core nucleon single-particle energy 
levels are in correct order. But, whenever we get the lower single-particle energy level for the 
lsi/2 state than those for the lpi/2, IP3/2 and 2si/ 2 states, the core neutron lsi/2 energy level 
jumped to become equal or larger than that of the 2si/ 2 , and we discarded the results in such 
cases. In this manner, the l«i/2 state in ^,0, | 9 Ca and |oZr could not be found. 

To get greater insight as to what happens, as an example, we show in Figs. 2 and 3 
the potential strengths and baryon density (together with the baryon density of nucleons), 
respectively for S° in |oZr, giving incorrect-" lsi/2(?)" (notation used in the figures) and correct- 
2si/2 states. "Pa.+CC." stands for potential for the effective Pauli blocking plus the channel 
coupling effects. By the terminology, incorrect-" lsi/ 2 (?)" state, we mean by that the obtained 
solution is incorrect, since both the single-particle energy for the E° and the core nucleon energy 
levels are incorrect. In this case, the obtained single-particle energies for the S° are, -7.4 and 
-8.1 MeV, respectively for the "lsi/2(?)" and 2si/2 states. Looking at the potential strengths 
in Fig. 2, the difference between the two states are only slight, and nearly indistinguishable. 
However, looking at the E° baryon density (scaled by a factor 5) in Fig. 3, it is clear that both 
solutions have two nodes, thus actually they are 2si/ 2 states. The S° density distributions are 
very different, because the core nucleon density distributions (particularly around 3 fm from the 
center of the hypernucleus) are appreciably different due to the incorrect core nucleon single- 
particle energy levels. The main difference in the calculation of the S c hypernuclei from that of 
the £ hypernuclei is the bare (effective) mass, where mass of the S c is much larger, and the wave 
function (if obtained) for the lsi/2 is expected to be localized considerably in the central region 
of each E c hypernucleus when the isospin dependent strong interaction and Coulomb forces are 
neglected. This may also be a consequence of an interplay between the isospin dependent force 
for the £° (the isospin third component is —1), and baryonic density in the central region. To 
test this we have repeated the calculation for |oZr single-particle energies by switching off both 

the Pauli blocking and the channel coupling effects simultaneously by setting As c = 0. The 
results are shown inside the brackets in Table 5. First, the energy levels from lp3/2 to 1^9/2 
states seem to be correct. But again the correct lsi/2 state solution is not found, in the sense 
that the lowest single-particle energy for s-state obtained is -23.0 MeV, and the corresponding 
wave function has two nodes (2si/2 state). Then, we have looked into the behavior of p meson 
field, particularly near the center of y Zr, which may show some hints. Indeed, we have found 
following, different and distinct behaviors of the p meson field between the 2si/ 2 state and that 
for the other (non-s) states: 

• For all the non-s states calculated, p meson field is repulsive and varies in the range 
+2 ~ +5 MeV up to ~ 4 fm from the center of |oZr, and gradually dies out towards the 
surface. 
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• For the 2si/ 2 state, it is strong attraction relative to that for the non-s states, —28 ~ —20 
MeV up to ~ 0.4 fm from the center of |oZr, and changes to repulsion around ~ 1 fm, to 
become maximum of ~ +7 MeV around 4 fm from the center, and similarly for the other 
states, gradually dies out towards the surface. 

This implies that the lsi/2 state would need much stronger attractive p meson field than 
that for the 2si/ 2 state to exist near the center of |oZr. However, such a case would certainly 
disturb the core neutrons which are bound deeply nearer to the center of |oZr, and result in 
giving incorrect core neutron energy levels as we have to discard such a solution observed. 
Therefore we discard such a solution. Thus, we would conclude that the peculiar feature of not 
finding the correct lsi/2 state in |oZr (and probably also for | 9 Ca) is due to the behavior of 
the isospin-dependent p meson field at very near the center of |oZr. If one ignores the incorrect 
energy levels of the core neutrons, one would get a solution and regard it as a "lsi/2" state, 
which is not correct in reality. As we will discuss later again in case of lsi/2 state in |° 9 Pb, 
indeed the isospin dependent p-meson field is expected to work in a peculiar manner for the 
E°-lsi/2 state in neutron rich TP C hypernuclei. 

On the other hand, recent experiment [49] suggests that a strongly repulsive E-nucleus 
potential with a non-zero size of the imaginary part is favorable to reproduce the measured 
spectra for middle and large baryon number hypernuclei. If we wish to be consistent with this 
result, we need to introduce even stronger repulsive potential, either by the channel coupling 
effect, or based on the quark model analysis to reproduce the E-atom data without resorting 
to the channel coupling effect [50], for the treatments of both the E and E c hypernuclei. In 
this case, the repulsion should be strong enough so that it entails no bound state for the 
E hypernuclei (and probably also for the E c hypernuclei) in the middle and large baryon 
number hypernuclei. This will need a more elaborate investigation in the future, using further 
accumulated experimental data and analyses. 

Next, in Figs. 4 and 5 we show the potential strengths for the lsi/2 state found in ^Ca 
and ^ 09 Pb (j = E 9 ' + ,H 9 ). Recall that the effects are different for the A+ and E c as they are 
true for the A and E. As in the limit of nuclear matter in Section 2.2, the scalar and vector 
potentials for the heavy baryons are also quite similar to the strange hyperons with the same 
light quark numbers in the corresponding hypernuclei. This feature also holds between the 
E c and E in the corresponding hypernuclei, except for a contributions due to the differences 
in charges, where the Coulomb force affects the baryon density distributions and then the 
vector and scalar potentials are also slightly modified. Thus, as far as the total baryon density 
distributions and the scalar and vector potentials are concerned, among/between the members 
within each multiplet, (A, A+,A&), (E, E c ) and (S, S c ) hypernuclei, they show quite similar 
features. However, in realistic nuclei, the Coulomb force plays a crucial role as mentioned 
before, and the single-particle energies thus obtained show very different features within each 
hypernuclei multiplet. Of course, the mass differences within the multiplet is also the dominant 
source for the differences in the obtained single-particle energies. 

The probability density distributions for the heavy baryons for the lsi/2 state, in ^Ca 
and ^ 09 Pb (j = E° ,+ ,H 9 ), are shown in Fig. 6. (We mention that the wave function for the 
lsi/2 state for |oCa is correct, because it has one node as shown in Fig. 6.) It is interesting to 
compare the E 9 and S 9 probability density distributions. Due to the Pauli blocking and the 
channel coupling effects, the probability density distributions for the E 9 are pushed away from 
the origin compared to those for the S 9 in both |oCa and | 99 Pb. In particular, the E 9 density 
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distributions in |° 9 Pb are really pushed away from the central region, and thus nearly losing 
the character of a typical lsi/2 state wave function. This fact shows that the isospin dependent 
p-meson mean field for the E° lsi/ 2 state in central region of the neutron rich nuclei is very 
important as discussed before, and that explains why the correct lsi/2 state wave functions are 
not obtained in f? Ca and fjoZr, where the (isospin asymmetric) baryon density distributions in 
central region of these nuclei are expected to be larger than those for the |° 9 Pb as in the case of 
normal nuclei. (In the case of goO, we expect that the size of the nucleus is much smaller and 
E° feels more sensitively the isospin dependent p-meson mean field, due to the higher baryon 
density and the limitation of the mean field approximation than those for the |°, 9 Pb.) On the 
contrary, it is rather surprising that the S° probability density distributions are higher and 
localized in the central region than those for the S°' + , although one can naively expect that 
they have an opposite characteristic, because of the smaller scalar attractive potential for the 
Thus, in the present calculation the effects of the Pauli blocking, and particularly the 
channel coupling, play an important role for these different features. Then, we must be careful 
in drawing any definite conclusions from the E° and E+ probability density distributions. 

Finally, we show in Table 6 the calculated binding energy per baryon, —E/A, r.m.s charge 
radius, r c h, and r.m.s radii of the heavy baryon (also A), neutron and proton distributions, rj, 
r n and r p , respectively. The results listed in Table 6 are calculated with the lsi/2 heavy baryon 
configuration for all cases. At a first glance, it is very clear that the r.m.s radius for the Aj, is 
very small compared to those for other baryons within the same baryon number hypernuclei as 
one would expect, due to heavy mass. On the other hand, r c h and/or r p for the A& hypernuclei 
are usually the largest in the all hypernuclei calculated. This implies that the protons in the 
core nucleus are relatively more pushed away compared to the other hypernuclei, although such 
feature is not seen for r n . The radii, r ch , r n and r p , may be grouped in similar magnitudes 
for all heavy baryon hypernuclei and A hypernuclei within the same baryon number multiplet 
(^Ca and | 9 Ca in Table 6 may be grouped together), reflecting the fact that the effect of 
the embedded baryon on these quantities are of order, ~ M c ,a/[(A — 1)M N + M c ,a\- As for 
the binding energy per baryon, the energy of Af, hypernuclei is usually the largest among the 
same baryon number hypernuclei. One of the largest contributions for this is the single-particle 
energy of the lsi/2 state, even after dividing by the total baryon number. 

5 Summary and discussion 

In summary, we have completed systematic studies of the A+, E c , S c and A& hypernuclei in the 
QMC model, extending the study made for strange hypernuclei [9]. The spin-orbit potentials 
for the heavy baryon hypernuclei are negligible for the single-particle energies, because of their 
heavy masses. Our results suggest that the formation of A+, S° and A& hypernuclei is expected 
to be quite likely, while that of the E+ + and is unlikely. For the E° and E+ hypernuclei, 
although there may be some possibilities to be formed, it is difficult to draw a solid conclusion 
in view of the uncertainties in the model and approximations made in the calculations. 

For the E c and S c hypernuclei, the Coulomb force is crucial in the possibility to form the 
hypernuclei, in particular for the S+ ,++ hypernuclei. The results for the E c hypernuclei imply 
that the phenomenologically introduced channel coupling effects may be overestimated when 
using the same strength as that for the E hypernuclei. Thus, we need to investigate further how 
such effects can be included consistently with the underlying quark structure. In this sense, 
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Table 6: Binding energy per baryon, —E/A (in MeV), r.m.s charge radius, r c h, and r.m.s radii 
of the A and heavy baryons, rj, neutron, r n , and proton, r p (in fm) for j 7 0, ^Ca, | 9 Ca, 91 Zr 
and ^ 09 Pb (j = A, A+, E c , S c , A 6 ). The configurations of the A and heavy baryon j, are lsi/ 2 
for all hypernuclei. 



hypernuclei —E/A r c h r 



\ 7 6.37 2.84 2.49 2.59 2.72 

^ 7 + 6.42 2.85 2.19 2.58 2.73 

]?+0 6.10 2.83 2.57 2.60 2.71 

||0 6.01 2.81 2.34 2.61 2.69 

5.84 2.81 2.70 2.61 2.69 

Ifo 6.69 2.87 1.81 2.57 2.75 

fCa 7T58 3.51 2.81 3.31 3.42 

f+Ca 7.58 3.51 2.66 3.31 3.42 

ffoCa 7.54 3.51 2.79 3.31 3.42 

£+Ca 7.42 3.50 3.21 3.31 3.41 

||Ca 7.43 3.49 2.76 3.32 3.40 

£Ca 7.32 3.49 3.11 3.32 3.39 

^Ca 7.72 3.53 2.22 3.30 3.43 

"YCa 7l8 3.54 2.84 3.63 3.45 

f+Ca 7.54 3.54 2.67 3.63 3.45 

4 $+Ca 7.39 3.54 3.17 3.64 3.44 

^++Ca 6.29 3.57 3.62 3.71 3.47 

| 9 oCa 7.34 3.53 2.57 3.65 3.43 

^Ca 7.32 3.53 3.19 3.65 3.43 

fCa 7.64 3.56 2.16 3.63 3.46 

~~ fZr 735 4.29 3.25 4.29 4.21 

91 +Zr 7.85 4.29 3.34 4.30 4.21 

91 +Zr 7.77 4.29 4.03 4.30 4.21 

||Zr 7.78 4.28 3.01 4.30 4.20 

fZr 7.93 4.30 2.63 4.29 4.22 

'f J Pb 7T35 5.49 3.99 5.67 5.43 

2° + 9 Pb 7.26 5.49 4.74 5.68 5.43 

||Pb 6.95 5.51 4.88 5.70 5.45 

|| 9 Pb 7.24 5.49 4.44 5.68 5.43 

|° 9 Pb 7.48 5.49 3.45 5.66 5.43 
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some of the results for the S c hypernuclei need to be taken with caution. 

Although there can be numerous speculations on the implications of the present results 
we would like to emphasize that our calculations indicate that the A+,S°, and A b hypernuclei 
would exist in realistic experimental conditions, but there may be lesser possibilities for the S° 
and hypernuclei. Furthermore, it is very unlikely that the £+ + and hypernuclei will be 
formed. Experiments at facilities like JHF would provide quantitative input to gain a better 
understanding of the interactions of heavy baryons with nuclear matter. Experiments at collid- 
ers such as RHIC, LHC and Fermilab could provide additional data to establish the formation 
and decay of such heavy baryon hypernuclei. A combination of these data inputs and a careful 
analysis, with the present calculations being considered as a first step, would give a valuable 
information about the physical implications for the presence of heavy quarks in finite nuclei or 
dense nuclear matter. Furthermore, the experiments would indicate certainly the shortcomings 
of the QMC model and in addition would provide a guide to a proper and consistent model to 
treat baryons with heavy quarks in finite nuclei or nuclear matter. Eventually we hope that 
these results will have important implications for studies in astrophysics of neutron stars and 
nuclear matter at very high densities. 

As for the other aspects of interests for the properties of heavy baryons in nuclear medium, 
additional studies are needed to investigate the semi-leptonic weak decay of A+, E c , S c and A b 
hyperons in nuclear medium. This will provide information on the weak current form factors 
of heavy baryons in nuclear medium, whether or not the axial coupling constant of such heavy 
system will be modified appreciably, as is the case implied for the quenching qa of nucleon, 
possibly due to the internal structure change of the heavy baryon in medium [51]. To determine 
the role of Pauli blocking and density in influencing the decay rates as compared to the free 
heavy baryons would be highly desirable. Such studies can provide important information on 
the hadronization of the quark-gluon plasma and the transport of hadrons in nuclear matter at 
high density. Will the high density lead to a slower decay and a higher probability to survive 
its passage through the material ? At present the study of the presence of heavy baryons in 
finite nuclei and/or nuclear matter is in its infancy. Careful investigations, both theoretical 
and experimental, would lead to a much better understanding of the properties of heavy quarks 
and/or heavy baryons in finite nuclei and nuclear matter. 
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Figure 1: Scalar potentials for hyperons, low-lying charmed and bottom baryons in symmetric 
nuclear matter taken from Ref. [2]. 
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Potential strengths for L c in E ° c Zr 

150 




Figure 2: Potential strengths for S° ,+ in |oZr. "Pa.+CC." stands for potential for the effective 
Pauli blocking plus the channel coupling effects. See the text for the meaning of "lsi/ 2 (?)". 
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Baryon densities in z ° c Zr 
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Potential strengths in c Ca 
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Figure 4: Potential strengths for the E°' + and S° in ^Ca (J = H°) for the lsi/ 2 state. See 
also caption of Fig. 2. 
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Potential strengths in c Pb 
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Figure 5: Potential strengths for the E° and S° in | 09 Pb (J = E°, H°) for the lsi/2 state. See 
also caption of Fig. 2. 
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Heavy baryon densities in c Ca and c Pb 
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Figure 6: and S° baryon (probability) density distributions for the \s\/2 state in ^Ca and 

209p b(j = E 0,+ £0). 
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